Pairing Fluctuations in The Attractive Hubbard Model in the Atomic Limit 
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BCS theory accounts for the pairing instabihty in the weak coupling limit, but fails to describe 
pairing fluctuations above Tc. One possibility for describing these fluctuations in the dilute limit 
is the T-matrix approximation. We critically examine various degrees of self-consistency in the T- 
matrix formalism, along with a non-diagrammatic two-particle self-consistent (TPSC) formulation, 
in the strong coupling regime, where an exact solution is readily available. We find that one 
particular degree of self-consistency is quite accurate, particularly at low temperature as evidenced 
by examining both static and dynamic properties. 

PACS numbers: 



In a recent essay titled, "Brainwashed by Feynman" 
, P. W. Anderson wrote about the different approaches 
to solving many-body problems. His main point was that 
the solution to some problems requires a leap beyond 
the "tried and true procedures" i.e. perturbation theory. 
He even argues that attempts to go beyond perturbation 
theory, for example by "summing all the diagrams" , are 
more or less futile. We generally agree with these senti- 
ments. 

On the other hand, we would argue when such a 
creative step occurs, as happened when BCS formu- 
lated their famous wavefunction for the superconducting 
ground state , there nonetheless remains a desire to re- 
cast the solution in terms of diagrams. In this particular 
case precisely this was accomplished shortly thereafter, 
by Gor'kov 3] and others. That this was indeed useful 
is clear by the subsequent work that followed concerning 
generalizations of BCS theory, such as Eliashberg the- 
ory ^j. The work by Thouless Q elucidated the nature 
of the superconducting instability in terms of an infinite 
set of diagrams — the ladder sum in the particle-particle 
(Cooper) channel. This allowed one to study the possibil- 
ity of a superconducting transition away from the weak 
coupling limit for which BCS was initially designed, and 
this indeed has occurred over the decades since, as ev- 
idenced by, for example, Refs. [IJl H H E|- Many 
more references are cited in Ref. All of this work 

is characterized by an attempt to cast the superconduct- 
ing instability in terms of a subset of Feynman diagrams; 
they thus suffer from the "Brainwashed by Feynman" 
(BBF) attitude documented in Ref. 1]. 

We should note that in recent years some work at- 
tempts to tackle the problem in a way which can- 
not be simply recast in terms of Feynman diagrams 
[l^ IT^ 113. Il5| . This approach has been reasonably 
successful as jud ged by comparison with (exact) Monte 



judged by c 

iialilEj. 



Carlo results PJ,! 

The purpose of this paper is to examine this problem 
once again, in a limit for which an exact solution is read- 
ily available, i.e. the atomic limit. This limit allows us 



to discriminate between candidate theories; we find that 
one in particular emerges as very accurate, particularly at 
low temperature. For simplicity we will use the attractive 
Hubbard model (AHM). The atomic limit is essentially 
in the strong coupling limit of this model; this is where 
the BBF approach should encounter the most difficulties. 
Nonetheless we bravely forge ahead, using the T-matrix 
approach pioneered by Thouless and others In the 
BBF approach one of the questions that researchers have 
been addressing is: what level of self-consistency (if any) 
most properly reproduces the correct result for the super- 
conducting instability over all regimes (weak to strong 
coupling) ? 

In what follows we will formulate t he p roblem in a 
manner which avoids Hartree diagrams hj. This mini- 
mizes the actual number of diagrams that need be con- 
sidered, and avoids the plethora of possible degrees of 
self-consistency. We also consider the two-particle self- 
consistent approach of Vilk and Tremblay (VT) ) as 
implemented by Kyung et al. for the AHM. The 

VT approach apparently works best for intermediate cou- 
pling. Remarkably it and one of the self-consistent T- 
matrix formulations both become exact at zero temper- 
ature in this strong coupling limit. 

Let us begin with the Hamiltonian. It is written |l7j 



H 



-|C^lE("'^T-«/2)Ki-«/2), (1) 
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where the first term is the usual hopping term for an elec- 
tron, the second is the chemical potential term, and the 
third describes the attractive interaction of strength \U\ 
between electrons on the same site. Note that the mean 
field expectation value of the electron density of a given 
spin species in the paramagnetic state, < n^o- >= no- — 
n/2, is subtracted from the electron density operator in 
the interaction term. Similarly, the modified chemical 
potential, /i', is given in terms of the actual chemical 
potential, fi, by the relation ji' = i.i + n\U\/2. The use 
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of this Hamiltonian allows us to ignore all Hartree dia- 
grams; they have been included automatically by using 
IjI instead of /i. 

The exact solution to this Hamiltonian is easy to ob- 
tain in the atomic limit (i.e. |i7| >> t , so that the prob- 
lem reduces to a single site problem) [l^: We find 



Gt(z) = 



1 - "i 



z + ii'~\U\n^ z + M' + |C/|(l-nx)' 



(2) 



and therefore, with the non-interacting limit, Go{z) = 
l/(z + ^'); ^ reference state, 



\U?n^{l-n{) 
z + /i' + |C/|(l-2nx)- 



(3) 



Note the simplicity of the exact solution. Accounting 
for the Hartree term in the modified chemical potential, 
the single electron propagator contains two poles, analo- 
gous to a lower and upper Hubbard band in the repulsive 
model, with energies separated by \U\. 

The electron density n is easy to obtain in terms of the 
modified chemical potential /i'. Inverting this relation, 
one obtains 



(l — n) — Tin (a(n)/n) 



(4) 
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FIG. 1; (a) Exact solutions and (b) (GoGo)Go results for ^! 
vs. temperature for different densities. In (b) the shaded 
region is the so-called Thouless region. Note that it is absent 
in the exact solution. 



Here the zero temperature result is given in the first term, 
and is quite simple, and finite temperature corrections 
are given by the second term, with 



a{n) = (1 - n)e-'^l^l/2 + sj n{2 - n) + {I - nfe-P\U\ . 

(5) 

A side benefit of using ^' is that the solutions are spread 
out as a function of electron density, which would not be 
the case if /i were used instead. Some of these solutions 
are shown in Fig. (fa). We also examine the two particle 
correlations. One may derive 

< n^n^ >= (f - n)N[-2^I' - \U\{1 - n)] , (6) 

where N{x) = l/{exp{(3x) — 1) is the Bose function. 
With the result of Eq. Q this 'simplifies' to 



2 n + b{n)exp{-/3\U\/2)'' 



(7) 



and&(n) = (f-n)e-'^l^l/2 + ^n(2 ~ n) + {1 - nfe-f^\^\. 
Again, as T ^ the result is quite simple and expected, 
< n^^ny >^ n/2. 

The T-matrix approximation to this problem results in 
a self energy of the form 0, 0, 0| 
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\U\xoii'^e) 



(8) 

where = inT(2m — 1) and ivi = i2'!TTi are the 

Fermion and Boson Matsubara frequencies, respectively, 



and m and i are integers. The 'bare' susceptibility 
Xo{ii^7i) is given by 

Xo(«i^n) = 4 X! Ga{iuJ„i)Gb{~iUJm + iVn) , (9) 
m 

where the subscripts a, 6, and c in the above two equa- 
tions can either be absent (to indicate that the fully in- 
teracting Green function should be used) or can be set 
to '0', to indicate that the non- interacting Green func- 
tion is used. In the former case, one is required to use 
Dyson's equation, G(iLOrn) = ^/{i^m + A*' ~ S(*'^m)) to 
self-consistently determine G and E. While the literature 
has some discussion of these various approximations, our 
purpose here is to examine them critically in the atomic 
limit. 

The non-self-consistent (NSC) approximation was first 
examined in two dimensions for the AHM in Rcfs. 
[iiil l2ll|. Fig. (lb) shows the results in the atomic limit. 
Comparison with Fig. (la) shows that this approxima- 
tion is poor: indeed the same conclusion occurs in two 
dimensions A 'Thouless region' in the fj! —T plane 
occurs; this is where an instability would occur if fi' was 
held fixed while the temperature is lowered; it is indicated 
in Fig. (lb) by the shaded region. The more physical 
procedure is to keep the electron density fixed. Then /i' 
adjusts so that a finite temperature transition is avoided. 
As in higher dimensions, the electrons form bound states, 
the Fermi sea is depleted, and a finite temperature tran- 
sition is avoided for the wrong reason^ as inspection of 
Fig. (fa) shows. 
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FIG. 2: A comparison of different approximations at n = 0.3 
using (a) and (b) energy. Note that the VT result (first 
iteration only) and the (GGo)Go T-matrix approximation are 
the most accurate; they both become exact as T 0. The 
second iteration of the VT result (not shown) is poorer than 
the first. The insets focus on the TPSC and (GGo)Go results 
at low temperatures; note in (b) how well the (GGo)Go curve 
approaches the exact result over an extended temperature 
range. 



We have examined all the possible variations for 
Eqs.lIHl ini) at low electron densities. The results for /x' 
vs. r at a low electron density are shown in Fig. (2a). 
In addition we have also computed the result for the 
Vilk-Tremblay ^ theory. At very low electron densities 
all approximations work reasonably well. At low (but 
not too low) electron densities, two emerge as particu- 
larly accurate, the {GGo)Go T-matrix theory, and the 
Vilk-Tremblay theory. We use a nomenclature for the 
T-matrix theory to correspond to {GaGb)Gc, where the 
a, b, and c refer to the labels in Eqs.ljHUni- 

Clearly, the Vilk-Tremblay result is most accurate. 
However, particularly at low temperatures, it is evident 
that the [GGq)Gq T-matrix approximation is also quite 
accurate. In fact, at lower electron densities (not shown) 
the (GGo)Gq result becomes more accurate than the VT 
result, and these two remain the front-runners for all 
electron densities. Interestingly, both the VT and the 
(GGo)Go results appear to become exact at low temper- 
atures. Indeed, in the {GGq)Go formulation the ansatz 



A2 



(10) 



emerges if one considers only the term in Eq. ((SJ with 
£ = 0. Then = T\U\^xoiO)/i^ - |C/|Xo(0)). In- 



sertion of this self energy into the one electron Green 
function allows one to evaluate the bare susceptibility 
in Eq. at zero frequency. Earlier work [ijj indi- 
cated that any degree of self-consistency in the 'bare' 
susceptibility drives the superconducting transition (sig- 
nalled by 1 = |i7|xo(0)) to zero temperature. Adopt- 
ing this requirement in this case gives the parameter 
in terms of \U\ and ^' . The number equation, 
n = ^^^G{ioJm) exp (iuJmO^) provides a second rela- 
tion between A^ and fi'; with these two equations all 
the zero temperature properties can be obtained analyt- 
ically, and, remarkably, they agree with the exact solu- 
tion. Thus, the suppression of the Thouless instability 
to zero temperature is a feature in common with higher 
dimensional solutions. For the VT theory, these results 
are for the first iteration . If one presses further with 
the second iteration, the results deteriorate, consistent 
with the observation of these authors . 

We can press further these comparisons, and in partic- 
ular probe two-particle correlations. An easy way to do 
this is by examining the energy per lattice site E, which 
normally would be simply proportional to the double oc- 
cupancy in the atomic limit. However, because of the 
special form of the Hamiltonian, Eq. the precise 

relation is 



E 



-\U\ < n^n^ > +\U\{ 



n^2 



We use the exact relation 

£; = i ^ T.{iujjn)G{iujm) exp (ia;™0+), 



(11) 



(12) 



to obtain the energy at any temperature, in the T-matrix 
theories. Fig. (2b) shows the energy determined within 
the various approximation schemes along with the exact 
and the VT result, for n = 0.3. Once again the VT and 
(GGo)Go results are exact at zero temperature. However, 
the VT result deviates immediately for T > 0, while the 
(GGo)Go result follows closely over some temperature 
range. That this latter result is true follows also from 
the analytical work described above. 

This agreement at T (note that neither approxi- 
mation is particularly accurate at intermediate temper- 
atures) raises the question of the origin of the effective 
potential, |C/pp|, in the VT theory. Within a conserving 
approxmiation, a frequency-independent irreducible ver- 
tex, I Upp I , is accompanied by single electron propagators 
that include a simple Hartree-like renormalization |l4l| . 
The (GGo)Go T-matrix approach, however, suggests that 
the origin of an effective interaction is from the (partial) 
renormalization of the single electron propagator. That 
is, we can construct a theory that appears similar to the 
VT theory, i.e. with unrenormalized propagators every- 
where in Eos. 1)819(1 . but with an effective interaction ver- 
tex. Eq.© suggests this will be accomplished by 



\U\xaiii^n) = lUefflxooiii'n), 



(13) 
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FIG. 3: The spectral function for the {GGo)Go T-matrix ap- 
proximation (symbols) , compared with the exact result (solid 
curve). The remarkable agreement shows that only partial 
self-consistency is a requirement to reproduce the exact re- 
sult. The VT result (both first and second iterations) are 
poorer in comparison. We have used a small artificial broad- 
ening, S = 0.01|L''| to plot these results. 



particle self-consistent (VT) phenomenology. We found 
that minimal self-consistency (the (GGo)Go theory pro- 
moted in particular in Ref. [l^) and the VT calculation 
both reproduced the exact result best, particularly at 
low temperatures, where both become exact. Surpris- 
ingly, 'improving' the degree of self-consistency within 
a T-matrix approach leads to less accurate results p^. 
Two particle correlations (summarized in the total energy 
in this work) are faithfully reproduced by the (GGo)Go 
calculation for a range of low temperatures, and, finally, 
the one electron spectral function is remarkably accu- 
rate, even at nonzero temperature. These results suggest 
that, at least for low electron densities, this T-matrix for- 
mulation is useful for higher dimension calculations. It 
would also be interesting to see how the various vertex 
corrections provide partial cancellation to the omitted 
self-consistent contributions. 



where Xooiivn) = 1 1]„ Go(iwm)Go(-iwm + ii^n)- Note 
that this requires |C4//| to depend on Matsubara fre- 
quency, but, in the spirit of VT (and many Parquet-like 
treatments of higher order corrections |23l|). we will use 
Eq. H13|l at zero Matsubara frequency to define an effec- 
tive potential. Comparisons with \Upp\ as obtained in the 
VT formalism show quantitative discrepancies, particu- 
larly as the temperature approaches zero. Thus, the two 
theories differ more substantively than Eq. (fT!^ would 
suggest. 

To probe further the analytical structure of the var- 
ious approximations, we examine the spectral function 
{A{uj) = — ImG(cj + id) /it). As already noted, the 
(GGo)Go reproduces the exact result at zero tempera- 
ture. This is shown in Fig. 3, where the spectral function 
for the (GGo)Go has been obtained through Pade ap- 
proximants. The agreement with the exact result at low 
temperatures is remarkable. On the other hand, the VT 
result in the first iteration gives a single pole, which is 
clearly inaccurate. The result obtained from the second 
iteration does contain two poles. However, as shown, it 
is fairly inaccurate. Thus, even though both the VT the- 
ory and the (GGo)Go T-matrix approximation give exact 
results for the chemical potential and the energy at low 
temperatures, only the latter fully reproduces the exact 
result as a function of frequency. Inspection of Fig. 2 
shows that 'improving' the degree of self-consistency de- 
teriorates the agreement. This would mean that at low 
temperatures some cancellation occurs between the fully 
self-consistent T-matrix diagrams ((GG)G theory) and 
the omitted vertex corrections. This possibility remains 
to be shown, however. 

In summary we have critically examined various ap- 
proximations for the attractive Hubbard model in the di- 
lute and strong coupling limit. We have used a variety of 
self-consistent T-niatrix approaches along with the two- 
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